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1. Introduction 

The Svetitsky-Yaffe conjecture [[|, ||] states that the Yang-Mills finite temperature transition in di- 
mension d + 1 is described by an effective spin model in d dimensions with short range interactions. 
Combining this idea with strong coupling expansions and inverse Monte-Carlo (IMC) methods we 
analyse the relationship between SU (3) YM theory in 3 + 1 dimensions and effective theories for- 
mulated as Z3 spin models in 3 dimensions. 



2. SU (3) and characters of representations 

Our effective operators are class functions on SU (3). With group elements in diagonal form, g = 
diag(e"^ ,e"^ ,e~ l ^ i+ ^), we associate a group character in the fundamental representation by 

& = *8 = Xw(g) = ^ + e < fe + e- i <*+« , (2.1) 



with the typical example being the Polyakov loop. The parameterisation ( |2.1| ) implies the reduced 
Haar measure on the maximal Abelian torus, 

dpL m d=J 2 d<$>\d<fo, J 2 = 15 -6X11 + 3^30 + 3^03-^22- (2.2) 

Using Young tableaux one can express all characters % pq with Dynkin labels [p,q] in terms of the 
fundamental ones, J 2 and . 



3. Observables 

We discuss YM theory on a N$ x A^-lattice. The Polyakov loop & x is measured in terms of its 
lattice average, 

' X 

The observable relevant for the analysis of antiferromagnetic phases is 

M^£^sgn(jt), sgn(jc) = (-l) E ^ (3.2) 

' X 

and measures the difference of the Polyakov loop on odd and even sublattices. 

Since we will have to deal with phases where the traced Polyakov loop is located halfway 
between the SU (3) center elements we project the value of the traced Polyakov loop onto the 
nearest 2^-axis and define a rotated Polyakov loop by 
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4. Effective models for Yang-Mills theory 

We start with the well-known lattice Wilson action 

and perform a strong coupling expansion (for small /3). Since the resulting 'operators' (Polyakov 
loop monomials) are dimensionless there is no natural ordering scheme. We therefore use a trun- 
cation scheme based on: 

• Ordering by powers of j8 which are closely related to the dimension of the corresponding 
group representations. 

• Ordering by the distance across which the Polyakov loops are coupled. 
In compact form the strong coupling expansion is given by 

Seff = I £ £ 4i (P) n S *iA = L ( 4 - 2 ) 
r Mi..M r l\..l r i=l i 

with the basic building blocks 

S^=xM^x)xU^y) + ^ (4-3) 

Here r counts the number of link operators contributing at each order. The coefficients c^j " ^ 
are the couplings between the operators S^.g, sitting at nearest-neighbor (NN) links l[ = (jc,-,y,-) 
in representation St\. The effective action hence describes a network of link operators that are 
collected into (possibly disconnected) 'polymers' contributing with 'weight' c^" ^ . One expects 
the 'weights' or couplings to decrease as the dimensions of the involved representations and inter- 
link distances increase. In a strong coupling (small j8) expansion truncated at £?(j3 ') one has r < k 

and the additional restriction h \M r \ <k with \&\ = p + q for a given representation 2% 

with Dynkin labels [p,q\. 

5. A toy model - mean field vs. Monte-Carlo 

We consider the St/ (3) model £| 

S = AiE(*ioW;ifoi(^)+cx.)+^ (5.1) 

(xy) {xy) 

A mean field approximation can be applied to approximately determine the associated phase 
diagram. We use the following ansatz for the distribution p of the field 

p[0\^Prt[&\=UP>&*) with P*^ = i M Z X l :Sgn{X \ = l ■ (5-2) 

[Poi&x) :sgn(jf) = -l 

The resulting phase diagram is displayed in Fig. |l| (left panel). 

A straightforward Monte-Carlo simulation on an 8 3 -lattice with a Metropolis algorithm using 
our jenLaTT package leads to a phase diagram (Fig. [I], right panel) similar to the one obtained 
by the mean field analysis. This agreement is due to the presence of a tri-critical point implying 
an upper critical dimension of three. In summary, the full phase structure consists of a symmetric 
phase (in the center of each panel of Fig. [j]), a ferromagnetic phase (upper left), an anti-center 
phase (lower left) and an antiferromagnetic (lower right) phase. The anti-center phase is related to 
the 'skewed' phase of [||]. 
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6. Algorithms for first and second order transitions 

The simulations for the microscopic YM theory were done using standard heat bath algorithms. For 
the effective models we employed standard Metropolis updates to find the phase diagram. In the 
vicinity of phase transitions we made use of the following specially designed update scheme. 

For first order transitions we used a multicanonical algorithm ^ improving the transition rate 
near critical points. For larger lattices the distribution p of the order parameter (denoted £) was 
predicted using the scaling relation 

logp(£,V)^A(£)+C(l)V . (6.1) 



For second order phase transitions algorithms of Wolff or Swendsen-Wang type lead to strong 
suppression of the dynamical critical exponent. These algorithms are useful for systems with invo- 
lutory global symmetries, where the local application of such symmetries leads to ergodic behavior 
of the system. For our system there is no ergodic symmetry. So we had to modify the well-known 
Wolff cluster algorithm [[7|] as follows: 

1. Choose a random number Nm between and V =N 3 . 

2. Do Nm standard Metropolis sweeps at randomly drawn lattice points. 

3. For a suitable fixed number N c \ repeat the steps for building a cluster by using the complex 
conjugation symmetry and its ^-symmetric equivalents. 

4. Do V — Nm additional Metropolis sweeps, again at randomly chosen lattice sites. 



7. Critical exponents for the antiferromagnetic SU (3) model 

For the model ( |5T| ) with A4 = we observe a second order transition between symmetric and 
antiferromagnetic phase. Critical exponents v and 7 may be introduced in terms of the relations 

c<N 1/v with U = l- ^ M \ , y=N 3 (M 2 ). (7.1) 



n \ wr/v du ( N M 
KM,c6t) oc N r ' , 



A Monte-Carlo simulation with our modified Wolff cluster algorithm leads to the following critical 
exponents in comparison to the Z3 Potts values: 
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exponent 


Z 3 Potts H 


minimal Polyakov 


V 


0.664(4) 


0.68(2) 


y/v 


1.973(9) 


1.96(2) 



As the exponents coincide (up to statistical errors) the SU (3) model is indeed in the same univer- 
sality class as the Z3 Potts model (the XY universality class). 



8. Inverse Monte-Carlo - the basics 

The inverse Monte-Carlo (IMC) method as designed in [Q] allows to determine (effective) actions 
from given configurations. In our case, these are Polyakov loops obtained from gauge configura- 
tions generated with the Wilson action. Via IMC we want to determine the couplings of truncated 
effective actions which (ideally) would give rise to the same distribution of Polyakov loop config- 
urations. 

The IMC procedure is based on an ansatz for the effective action of the type S e ff = L,-A,S ; . 
Translational invariance of the reduced Haar measure leads to Schwinger-Dyson equations (SDE), 
see below. They constitute an overdetermined linear system for the effective couplings A, which 
may be solved by least-square methods. As a further technical input we require a suitable normal- 
ization procedure to make sure that individual equations are appropriately weighted [|To|]. 



9. Geometric SDE from invariant group integrals 

Translational invariance of the Haar measure implies that 

JdHn aa r(g)(L a f)(g)=0 for/eL 2 (G) (9.1) 

with L a being the left derivative on the group. Choosing / = FL a % p with a class function F and a 
fundamental character % p [ 11 ] one obtains 

L a (FL" Xp ) = FL 2 Xp + (L a F)(L a Xp ), (9.2) 

and ( |9~I| ) reduces to 

= | ^Haar ^FL 2 Xp + ^aXp) . (9.3) 

Making use of L 2 %^ = —CnXfi an d of 

(LaZt)(L a Zv) = Ucn+c v )XnXv -\Y. C lvCxXx (9-4) 
z x 

with Casimir values and Clebsch-Gordan coefficients C^ v the equation can be specialized to the 
case of SU (3). For a suitable function F, the SDE finally become 

= /-y«,^ + (4^ z * - \^l)Si,^ + (6- \\^ z f)h^l) 

I 4 2 \ (9 " 5) 

-IM(4^; - + (6 - -\& z \ 2 )S L ^S m 
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10. Algebraic SDE 

For SU (3) (generalizations for SU (N) are possible [|TI[]) we have the identity 



/ d&d&* d&f = 



which holds for any function / vanishing on d£l. Choosing 

r)h 

/(^,^*)=/ 3 g(^,^*), 8x= " 
we obtain the 'algebraic SDE' 



, 3 dJ 2 z 
— 



(10.1) 



exp(— S), h = Si 



(10.2) 



(10.3) 



11. IMC results 

We have simulated the underlying YM the- 
ory with heat bath methods on different lat- 
tice sizes near the critical coupling. The cor- 
responding couplings for the effective mod- 



els were then determined via IMC [10]. The 




IMC codes were checked by simulating ef- 
fective theories with given input couplings 
which were consistently reproduced by the 
IMC procedure. In these tests the algebraic 
and geometric Ward identities led to compa- 
rable results, limited only by the statistical ac- 
curacy. 

We compared the Polyakov loop arising 
from simulations of full YM theory and from 
effective actions based on both geometric and 

algebraic SDE (on a 16 3 x4-lattice). We found that the algebraic identities outperformed the geo- 
metric ones in reproducing the YM critical behavior, in particular the critical coupling (Fig. ^). 

Simulations with algebraic SDE on a 16 3 x4-lattice allowed to determine up to 11 effective 
couplings as displayed in Fig. [3] (left panel). The dominant terms in the effective actions are 



8.0 



Figure 2: Results obtained from algebraic and geomet- 
ric SD equations compared to YM results. 



5i = £(Zlo(^koi(^y)+cc), S 3 



(xy) x 



(11.1) 



i.e. two NN hopping terms and one single-site ('potential') term. 

Finally we have extended the IMC procedure to deal with NN and next-to-NN terms up to order 
^fQ3 3Af t) m th e strong coupling expansion. This results in an unstable behavior in the rendering 
of observables which may be traced to the discontinuities associated with the first order phase 
transition (Fig. ||, right panel). 
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Figure 3: Couplings obtained for NN interactions up to G(p)^ Nt ) (left) and comparison of higher order 
(next-to-NN) effective theories (right). 



12. Conclusions 

SU (3) Polyakov loop models have a surprisingly rich phase structure when the effective couplings 
are allowed to vary unrestrictedly. Upon comparing critical exponents for the second-order antifer- 
romagnetic phase transition we have seen that the SU (3) Polyakov loop model is in the same uni- 
versality class as the Z3 Potts model. The near-perfect agreement between mean-field and Monte- 
Carlo results is due to the fact that the model has a d = 3 tricritical point. Matching the Polyakov 
loop models to SU (3) YM theory via IMC leads to stable results only for small lattices and a low 
number of couplings. Relaxing these restrictions leads to instabilities obscuring, in particular, the 
location of critical couplings. This behavior is due to the first order nature of the phase transition 



in lattice gluodynamics. Results for SU (2) YM theory [12], on the other hand, show that IMC is 
applicable for systems with second order transitions and leads to stable results. 
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